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THE EGNOG SURVEY: GAS EXCITATION IN NORMAL GALAXIES AT Z « 0.3 
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ABSTRACT 

As observations of molecular gas in galaxies are pushed to lower star formation rate galaxies at 
higher redshifts, it is becoming increasingly important to understand the conditions of the gas in 
these systems to properly infer their molecular gas content. The rotational transitions of the carbon 
monoxide (CO) molecule provide an excellent probe of the gas excitation conditions in these galaxies. 
In this paper we present the results from the gas excitation sample of the Evolution of molecular 
Gas in Normal Galaxies (EGNoG) survey at the Combined Array for Research in Millimeter-wave 
Astronomy (CARMA). This subset of the full EGNoG sample consists of four galaxies at z ^ 0.3 
with star formation rates of 40 — 65 yr~^ and stellar masses oi ^ 2x 10^^ M0. Using the 3 
mm and 1 mm bands at CARMA, we observe both the CO(J = 1^0) and CO(J = 3^2) 
transitions in these four galaxies in order to probe the excitation of the molecular gas. We report 
robust detections of both lines in three galaxies (and an upper limit on the fourth), with an average 



line ratio, r^i = ^co(3-2)/^co(i-o)' ^-^^ ^ ^•^'^ (with systematic errors < 40%), which implies 
sub-thermal excitation of the CO (J = 3 2) line. We conclude that the excitation of the gas in these 
massive, highly star- forming galaxies is consistent with normal star-forming galaxies such as local 
spirals, not starbursting systems like local ultra-luminous infrared galaxies. Since the EGNoG gas 
excitation sample galaxies are selected from the main sequence of star-forming galaxies, we suggest 
that this result is applicable to studies of main sequence galaxies at intermediate and high redshifts, 
supporting the assumptions made in studies that find molecular gas fractions in star forming galaxies 
at z ~ 1 — 2 to be an order of magnitude larger than what is observed locally. 
Subject headings: Galaxies:evolution — Galaxies: high-redshift — GalaxiesiISM 



1. INTRODUCTION 

In the past decade, molecular gas observations have 
begun probing the high redshift universe in a systematic 
way using increasingly powerful millimeter instruments. 
The picture that is emerging at redshifts 1 — 2 is similar 
in some respects to what we see in the local universe. 
Sub-mm galaxies (SMGs) are observed to be undergo- 
ing extreme starbursts as a result of major interactions 
or mergers (e.g. Sheth et al. 2004; Tacconi et al. 2008; 
Engel et al. 2010), equivalent to local ultra-luminous 
infrared galaxies (ULIRGs) (e.g. Sanders et al. 1986; 
Solomon et al. 1997; Downes & Solomon 1998). Star- 
forming galaxies at high redshifts akin to local spirals are 
becoming accessible as well. Recent works (Tacconi et al. 
2010; Daddi et al. 2010; see also work by Baker et al. 
2004 on a 2: = 2.7 Lyman break galaxy) suggest that 
z 1 — 2 star-forming galaxies (with star formation 
rates (SFRs) of 50 - 200 M© yr-^) are scaled up ver- 
sions of local spirals, forming stars in a steady mode (not 
triggered by interaction), despite hosting star formation 
rates at the level of typical local starburst systems like 
luminous infrared galaxies (LIRGs) and ULIRGS. 

While galaxies classified as LIRGs or ULIRGs (by their 
infrared luminosities only) have typically been associated 
with starbursting and merging systems by analogy to 
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galaxies in the local universe, it is becoming clear that 
this connection does not hold at high redshifts. Morpho- 
logical studies find that while 50% of local LIRGs show 
evidence for major mergers (Wang et al. 2006), that frac- 
tion appears to decrease toward high redshifts: Bell et al. 
(2005) find that more than half of intensely star-forming 
galaxies at z ^ 0.7 have spiral morphologies and fewer 
than 30% show evidence of strong interaction. Further, 
the typical SFR of (normal) star-forming galaxies (like 
local spirals) increases toward higher redshift. Star- 
forming galaxies have been observed to obey a tight re- 
lation between stellar mass and star formation rate out 
to z ~ 2.5 (e.g. Brinchmann et al. 2004; Noeske et al. 
2007; Elbaz et al. 2007; Daddi et al. 2007; Pannella et al. 
2009; Elbaz et al. 2011; Karim et al. 2011). This 'main 
sequence' defines what is normal for star- forming galaxies 
as a function of redshift, showing that typical SFRs in- 
crease with redshift. Throughout this work, we use SFG 
to refer to main sequence (normal) star- forming galaxies. 

The increase in the star formation rate is mirrored in 
the molecular gas fraction of these systems. While stud- 
ies of local spirals (e.g. FCRAO survey. Young et al. 
1995; BIMA SONG, Regan et al. 2001, Heifer et al. 
2003; HERACLES, Leroy et al. 2009; also Combes et al. 
1994; Kuno et al. 2007) find average molecular gas frac- 
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stellar mass), observations of high redshift SFGs suggest 
molecular gas fractions of 20-80%, an order of magni- 
tude higher than local spirals. It is clear that in order 
to understand these SFGs at z ~ 1 — 2, we must investi- 
gate the condition of the molecular gas which is forming 
stars at such an enhanced rate. Do these systems truly 
hold massive reservoirs of molecular gas or is enhanced 
excitation of the gas misleading the interpretation of the 
observations? 

The rotational transitions of the carbon monoxide 
(CO) molecule provide a direct probe of the excita- 
tion of the molecular gas in galaxies. Local starburst 
galaxies and ULIRGs (e.g. Bayet et al. 2004; WeiB et al. 
2005b; Papadopoulos et al. 2007; Greve et al. 2009) and 
high redshift SMGs and quasars (e.g. WeiB et al. 2005a; 
Riechers et al. 2006; WeiB et al. 2007; Riechers 2011; 
Riechers et al. 2011) show signatures of excited molec- 
ular gas: observed CO line spectral energy distributions 
peak at Jupper > 5 with thermalized lines up to Jupper > 
3. In contrast, studies of the Milky Way (Fixsen et al. 
1999) and local SFGs (e.g. Mauersberger et al. 1999; 
Yao et al. 2003; Mao et al. 2010) find a wide spread of 
excitation conditions, with an average that implies less- 
excited gas, where the Jupper = 3 line is already sub- 
thermal. This suggests that the CO (J = 3^2) line in 
particular is an indicator of the star formation character 
of a galaxy: 'starburst' versus 'normal'. 

More important still, the ratio of higher rotational lines 
of CO to the CO (J = 1^0) line is necessary to trans- 
late the observed CO line luminosity to a molecular gas 
mass using Xco or aco, since this conversion factor is 
calibrated for the CO (J = 1^0) luminosity. Therefore, 
the measurement of rji = ^co(j-(j-i))/^co(i-o) 
only informs our interpretation of the current interme- 
diate and high redshift CO line studies, but is critical 
in the era of ALMA, which provides an order of magni- 
tude increase in sensitivity, making the CO lines in high 
redshift galaxies more accessible. The CO (J = 3 2) 
line in particular, observed in the 1mm, 2mm and 3mm 
bands, probes the molecular gas in galaxies at z~ 0.3 — 3. 
This redshift range is of particular interest since it in- 
cludes the peak of the star formation rate density of the 
universe {z ^ 1 — 2) and therefore the height of galaxy 
building. 

While several studies have measured r^i in the lo- 
cal universe (Mauersberger et al. 1999; Yao et al. 2003; 
Mao et al. 2010; Papadopoulos et al. 2012), the measure- 
ment of line ratios in intermediate and high redshift 
galaxies has mostly been limited to SMGs and quasars 
(see the review by Solomon & Vanden Bout 2005). Our 
knowledge of CO line ratios in SFGs at 2; ~ 1 — 2 is lim- 
ited to one study at z = 1.5 (presented in two papers: 
Dannerbauer et al. 2009 and Aravena et al. 2010) which 
measures r2i in three galaxies and rsi in one galaxy. It 
is clear that more work is needed to better constrain the 
line ratios in intermediate and high redshift SFGs to in- 
terpret existing and future data and better understand 
the state of the molecular gas in these systems. 

As part of the Evolution of molecular Gas in Normal 
Galaxies (EGNoG) survey, we observe both the CO (J = 
1^0) and CO (J = 3^2) lines in four galaxies at 
z ^ 0.3 (the gas excitation sample), more than doubhng 
the number of SFGs at z > 0.1 in which CO line ratios 



have been measured. In this paper, we present the rsi 
values for the EGNoG gas excitation sample and compare 
to previous work at low and high redshift s. We discuss 
the implications of this measurement for the excitation 
of the molecular gas in the observed galaxies as well as 
for high-redshift SFGs. 

The paper is organized as follows: in Section 2 we give 
a brief description of the EGNoG survey as a whole, 
describe the selection of the gas excitation sample and 
present the redshifts, SFRs and stellar masses of the 
galaxies in this work; in Section 3 we describe the obser- 
vations and data reduction; in Section 4 we present the 
analysis of rsi in these galaxies; in Section 5 we discuss 
the implications of this work; and we give some con- 
cluding remarks in Section 6. The data reduction and 
measurement of fiuxes is discussed in detail in Appendix 
A and moment maps of the detected CO emission are 
presented in Appendix B. Throughout this work, we use 
a ACDM cosmology with (h, I^m, ^a) = (0.7, 0.3, 0.7). 

2. THE EGNOG SURVEY 

The EGNoG survey is a key project at the Com- 
bined Array for Research in Millimeter- wave Astronomy 
(CARMA).^ By observing rotational lines of the ^^CO 
molecule, the EGNoG survey traces the molecular gas 
in 31 galaxies from z = 0.05 to 0.5, where the sig- 
nificant evolution from gas-rich galaxies at z ~ 1 — 2 
(Daddi et al. 2010; Tacconi et al. 2010) to relatively gas- 
poor local galaxies (e.g. BIMA SONG, Heifer et al. 
2003; HERACLES, Leroy et al. 2009; CARMA STING, 
Rahman et al. 2012) remains almost entirely unobserved. 
The full EGNoG sample is split into 4 redshift bins: 
bin A, 13 sources, z = 0.05 — 0.1; bin B, 10 sources, 
z = 0.16 - 0.20; bin C, 4 sources, z = 0.28 - 0.32; and 
bin D, 4 sources, z = 0.47 — 0.53. The full survey is 
presented in the forthcoming paper Bauermeister et al. 
2013b, in preparation. 

In this paper, we present the gas excitation sample (bin 
C) of the EGNoG survey, the subset of the full survey 
sample for which we observe both the CO (J = 1^0) 
and CO (J = 3 2) lines. Observations were made 
using the CARMA 15-element array with the 3mm band 
(single-polarization) for the CO (J = 1^0) observations 
and the 1mm band (dual-polarization) for the CO (J = 
3^2) observations. 

2.1. Sample Selection 

The galaxies of EGNoG bin C are drawn from the main 
spectroscopic sample of the Sloan Digital Sky Survey 
(SDSS), Data Release 7 (York et al. 2000; Strauss et al. 
2002; Abazajian et al. 2009). The galaxies were se- 
lected from the parent sample to be as representative 
as possible of the main sequence of star-forming galax- 
ies. The main sequence is the tight correlation be- 
tween and SFR that has been observed over a large 
range of redshifts: e.g. z ^ (Brinchmann et al. 2004), 
z - 0.2 - 1 (Noeske et al. 2007), z - 1 - 2 (Elbaz et al. 
2007; Daddi et al. 2007; Pannella et al. 2009), (see also 
the summary of recent results in Dutton et al. (2010)). 

^ CARMA is a 3-band, 23-element millimeter interferometer 
jointly operated by the California Institute of Technology, Uni- 
versity of California Berkeley, University of Chicago, University of 
Illinois at Urbana-Champaign, and University of Maryland. 
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EGNoG 
Name 


SDSS identification 


RA 


Dec 


z 


log(M*/M0) 


SFR 
(M0 yr-i) 


CI" 
C2 
C3 
C4 


SDSS J092831. 94+252313.9 
SDSS J090636.69+162807.1 
SDSS J132047.13+160643.7 
SDSS J133849.18+403331.7 


09:28:31.941 
09:06:36.694 
13:20:47.139 
13:38:49.189 


+25:23:13.925 
+16:28:07.136 
+16:06:43.720 
+40:33:31.748 


0.283020 + 0.000022 
0.300622 ± 0.000010 
0.312361 ±0.000014 
0.285380 + 0.000015 


11 24+0-i° 
11 20+°-29 
11 46+°-^^ 
11 26+°-^^ 


Qo fT+SS.Q 
'^^•'-25.6 
p:+90.1 
219 

50.5ltt:l 



Table 1 

Basic information. Derived quantities (z, M*, SFR) are from the MPA-JHU group (see text). 
" indicates duplicate source in SDSS: the average value is reported for M* and SFR. 



Building on these observations, a few authors have 
attempted to describe the main sequence relation at 
< z < 2.5 with one equation (Bouche et al. 2010; 
Karim et al. 2011; Elbaz et al. 2011). In order to clas- 
sify the EGNoG galaxies, we adopt a relation which 
roughly agrees with the relations from Bouche et al. 

(2010) , Karim et al. (2011) and Elbaz et al. (2011) (see 
Bauermeister et al. 2013b, in preparation, for a complete 
description): 

/ M \-o.2 
sSFRusiGyr-') = 0.07(1 + zf-' (l^nij^j 

where sSFRms is the specific star formation rate (sSFR= 
SFR/M*) of the main sequence (MS) of star- forming 
galaxies. We differentiate 'starburst' (SB) from 'nor- 
mal' galaxies according the criteria of Rodighiero et al. 

(2011) : sSFRsB > 4 x sSFRus- Therefore, SFGs (nor- 
mal star-forming galaxies) lie roughly within a factor of 
4 of the main sequence sSFR (which is a function of 
and z) and starbursts lie at sSFR values larger than 4 
times the main sequence sSFR. 

In the selection of the EGNoG sample, we apply the 
following criteria in each redshift bin in order to identify 
non-interacting, star-forming galaxies lying as close to 
the main sequence as possible. Star- forming galaxies are 
selected using the BPT diagram (Baldwin et al. 1981) 
criteria from Kauffmann et al. (2003) (rejecting sources 
with active galactic nuclei). Obviously interacting galax- 
ies were excluded via visual inspection of the SDSS op- 
tical images. However, some interacting galaxies may 
be in the sample due to the difficulty of identification 
at the modest resolution of the SDSS images. Practical 
considerations imposed the following further constraints. 
We required a spectroscopic redshift so that the error 
in the redshift is small enough to ensure that CO emis- 
sion would be captured within the observed bandwidth. 
We excluded galaxies with SFRs below a minimum value 
estimated from the instrument sensitivity. 

The result of the sample selection is illustrated in Fig- 
ure 1, which shows stellar mass versus SFR of the EG- 
NoG bin C galaxies (red points). The blue shading indi- 
cates the (logarithm of the) density of points in the 
- SFR plane for all star-forming galaxies (spectroscopic 
targets only; ^ 1300 galaxies) at 2; = 0.25 — 0.35. We plot 
galaxies in a slightly larger redshift range than the bin 
C specification in order to better capture the behavior of 
the main sequence. The main sequence of star- forming 
galaxies at z = 0.3 is indicated by the solid black line, 
with the starburst cutoff indicated by the dashed black 
line. In this plot, the low-mass, low- SFR end of the main 
sequence is sparsely sampled due to the limited num- 
ber of spectroscopically targeted, low-SFR, star-forming 



galaxies in the SDSS at z ^ 0.3. While the EGNoG bin C 
galaxies are high-MH<, high-SFR galaxes, they lie within 
the expected scatter of the main sequence (as defined in 
Equation 1), and are not classified as starburst galax- 
ies according to the prescription from Rodighiero et al. 
(2011) (presented above). 

2.2. Redshifts, Stellar Masses and Star Formation Rates 

Spectroscopic redshifts are from David Schlegel's 
spZbest files produced by the Princeton- ID code, 
spe cBS^. The stellar masses and SFRs of galaxies 
in the SDSS DR7 are provided by the Max-Planck- 
Institute for Astrophysics - John Hopkins Univer- 
sity group (http://www.inpa-garching.inpg.de/SDSS). 
Stellar masses are derived by fitting SDSS ugriz photom- 
etry to a grid of models spanning a wide range of star for- 
mation histories. This method is found to compare quite 
well with the Kauffmann et al. (2003) methodology using 
spectral features (more detail on this comparison is found 
on the website above). Star formation rates are derived 
by fitting the fiuxes of no less than 5 emission lines using 
the method described in Brinchmann et al. (2004). Both 
stellar masses and star formation rates are derived using 
a Bayesian analysis, producing probability distributions 
of each quantity for each galaxy. These distributions for 




Figure 1. Stellar mass versus SFR in EGNoG bin C. Red points 
and error bars show EGNoG galaxies and blue shading shows the 
logarithm of the density (in the M*-SFR plane) of star-forming 
galaxies from the parent dataset at z = 0.25 — 0.35 (slightly larger 
than the EGNoG redshift range so that more points may be in- 
cluded to better capture the behavior of the main sequence). The 
solid black line indicates the main sequence (Equation 1) at the av- 
erage redshift of the bin. The starburst (SB) criterion is indicated 
by the black dashed line. 

^ See http://spectro.princeton.edu/ for more information 
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Figure 2. Probability distribution functions for the stellar masses 
(left panel) and star formation rates (right panel) from the MPA- 
JHU group for bin C sources. Data points give the 2.5, 16, 50, 
84 and 97.5 percentiles of the distributions. The median, mean 
and mode are indicated by the vertical solid, dashed and dotted 
lines respectively. In cases where a duplicate exists (source CI), 
the duplicate values are plotted in red. 

the four bin C galaxies are shown in Figure 2. We take 
the median of the distribution, with errors indicated by 
the 16th and 84th percentile points. In the case where a 
duplicate SDSS source exists (as a result of SDSS auto- 
mated source-finding) , we take the average of the two me- 
dian values and use the lowest (highest) 16th(84th) per- 
centile value to indicate the negative (positive) error. The 
redshifts, stellar masses and star formation rates (with 
errors) are given in Table 1. 

3. CARMA OBSERVATIONS 

The four bin C galaxies were observed in two rota- 
tional transitions of the CO molecule: CO (J = 1^0) 
Urest = 115.3 GHz, Uohs ^ 88 GHz and CO(J = 3^2) 
^rest = 345.9 GHz, i^obs ~ 266 GHz. At each fre- 
quency, each galaxy was observed over several different 
days. Each dataset includes observations of a nearby 
quasar for phase calibration (taken every 15-20 minutes), 
a bright quasar for passband calibration and either Mars 
or MWC349 for flux calibration (in most cases). 

The reduction of all observations for this survey was 
carried out within the EGN^ data reduction infras- 

^ http : //carma . astro . umd . edu/wiki/index . php/EGN 



tructure (based on the MIS pipeline; Pound & Teuben 
2012) using the Multichannel Image Reconstruction, Im- 
age Analysis and Display (MIRIAD; Sault et al. 2011) 
package for radio interferometer data reduction. Our 
data analysis also used the mi ri ad-python software 
package (Williams et al. 2012). The data were flagged, 
passband-calibrated and phase calibrated in the stan- 
dard way. Final images were created using invert 
with options=inosaic in order to properly handle and 
correct for the three different primary beam patterns. 
All observations are single-pointing. We describe the 
CO(J = 1^0) and CO(J = 3^2) observations indi- 
vidually below. A full description of the data reduction 
and flux measurement is given in Appendix A. 

3.1. CO(J = 1^0) 

The CO (J =1^0) transition for the galaxies pre- 
sented here lies in the 3 mm band of CARMA (single- 
polarization, linearly polarized feeds). The wide band- 
width of the CARMA correlator allowed us to observe 
both the i2C0(J = 1^0) and the ^^CO{J = 1^0) 
lines simultaneously. The ^^CO (^^CO) line was observed 
with five (three) overlapping 500 MHz bands, covering ^ 
7000 (5000) km s"^ total, at 42 km s"^ resolution, with 
3-bit sampling. These observations were carried out from 
August to November 2011 in CARMA's D configuration, 
with 11 — 150 m baselines yielding a typical synthesized 
beam of 4.8'^ x 3.9'^ at 88 GHz. Each galaxy was ob- 
served for 20 to 30 hours (time on-source) in moderate 
to good weather conditions for 3mm observation, yield- 
ing final images with rms noise of ~ 1.2 mJy beam~^ in 
a 42 km s~^ channel. The flux scale in each dataset is 
set by the flux of the phase calibrator, which is deter- 
mined from the flux calibrator. For these data, we use 
phase calibrators 0854+201 (4 Jy, 9% linearly polarized), 
1357+193 (0.8 Jy) and 1310+323 (1.7 Jy). 

All four galaxies were clearly detected in the CO(J = 
1^0) line, as shown in Figure 3. The top panels show 
the vector- averaged Real and Imaginary amplitudes and 
phase of the calibrated uv data versus velocity for the 
CO (J =1^0) line. For a compact source at the center 
of the field of view, the i?ea/ part shows the signal without 
a noise bias, and the Imaginary part shows the noise. In 
all four cases, we see coherent emission (larger Real am- 
plitudes coincident with noise-like Imaginary amplitudes 
and phases of ~ 0) over multiple velocity channels, in- 
dicative of a detection. We report no detection of the 
^^CO(J =1^0) line, which is expected to be weaker 
than the ^^00(7 =1^0) line by a factor of 7 to 17 
(Rickard & Blitz 1985). 

3.2. CO(J = 3^2) 

The CO (J = 3^2) transition at the redshift of 
the galaxies discussed here lies in the CARMA 1 mm 
band (dual-polarization, circularly polarized feeds). We 
again observed both the ^^CO(J = 3^2) and the 
^^CO(J = 3^2) lines simultaneously. The ^^CO 
(^^CO) line was observed with three (one) overlapping 
500 MHz bands, covering ^ 1500 (550) km s"^ total, at 
14 km s~^ resolution, with 3-bit sampling. These obser- 
vations were carried out in CARMA's E and D configu- 
ration during August 2011 and April 2012, respectively. 
Source C4 was observed entirely in E configuration, while 
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Figure 3. Vector- aver aged Real (solid) and Imaginary (dotted) amplitudes (Real,Imag in mJy beam ^) and phase (ph, degrees) of the 
calibrated uv data versus velocity (km s"-*^) for the CO(J = 1—^0) (red, top panel) and CO(J = 3—^2) (blue, bottom panel) transitions 
in each of the four sources (indicated in upper left corner). With the higher resolution of the CO (J = 3—^2) data, we observe the peak 
of the CO(J = 3—^2) emission to be offset (< 1.5^' ^ 6 kpc) from the centers of galaxies CI and C2. In these cases, we calculate the uv 
spectra at the peak of the CO(J = 3—^2) emission. For the CO(J = 1 — 0) line of galaxy CI, we average 2 channels together in order to 
increase the signal to noise ratio. Large Real amplitudes (the Imaginary part shows the noise) coincident with phases of ~ over multiple 
velocity channels indicate a detection. 



the other three sources were observed mostly in D con- 
figuration. The E configuration has 8 — 66 m baselines 
yielding a typical synthesized beam of 3.2" x 2.5" at 266 
GHz. The D configuration has 11 — 150 m baselines yield- 
ing a typical synthesized beam of 1.7" x 1.5" at 266 GHz. 
Each galaxy was observed for 2 to 6.5 hours (time on- 
source) in good weather conditions for 1mm observation, 
yielding final images with rms noise of 5-10 mJy beam~^ 
in a 42 km s~^ channel. For these data, we use phase 
calibrators 0854+201 (2 Jy in August 2011, 4 Jy in April 
2012), 1224+213 (0.6 Jy) and 1310+323 (0.6 Jy). 

Sources C2, C3 and C4 were detected at the ^ 5a 
level in the CO (J = 3 2) line. However, CI, with its 
wide velocity profile (observed in the CO (J = 1^0) 
line), was only marginally detected and we give only an 
upper limit on the CO (J = 3 2) fiux. The vector- 
averaged Real and Imaginary amplitudes and phase of 
the calibrated uv data versus velocity for the CO (J = 
3^2) line are shown in the bottom panels of Figure 3. 
The peak of the CO(J = 3^2) emission in galaxies CI 
and C2 is offset from the nominal center of each galaxy, 
so the uv spectra are calculated at this slightly offset 
(< 1.5'' ~ 6 kpc)) position. We report no detection in 
the ^^CO(J = 3 2) line for any galaxy. 

3.2.1. Source CI Upper Limit 

Due to the wide integrated velocity profile of source 
CI, the CO (J = 3 2) line was only marginally de- 
tected. The channel maps did not show evidence of a 
source upon visual inspection, but an integrated spec- 
trum made of a circular region 4.5'' in radius at the 
center of the image suggests a 3cr detection. The line 
fiux and other quantities are calculated from this spec- 
trum, over the velocities of the CO (J =1^0) line for 
this galaxy. These values should be taken as an upper 
limit on the true fiux. The error in the fiux measure- 



ment {Sco^ Jy km s~^) is calculated from the pixel noise, 
<jp (Jy beam""*^): crsco — y^^p/^^eq where SV is 

the channel width (km s~^), Np is the number of pixels 
summed and TVeq is the number of pixels equivalent to 
the beam. 



3.3. Derived Properties of the CO Emission 

Table 2 presents the quantities we derive from the CO 
emission images. The CO line luminosity is calculated 
from the hue fiux {Sco in Jy km s~^, calculated as de- 
scribed in Section A. 2) following 



L'co = 3.25 X 10^ ^co 



.-2 

obs 



.(1 + ^)" 



(2) 



(see the review by Solomon & Vanden Bout 2005), 
where Vohs is in GHz and rcom is the comoving distance 
in Mpc. The units of Lqq are K km s~^ pc^. We report 
the measurement error for Sco- The error reported for 
Lqq includes both the measurement error of ^co and a 
30% systematic error, added in quadrature (see Section 
A. 2 for more details). 

The center velocity, '^center, is the fiux- weighted aver- 
age velocity of the galaxy-integrated spectrum {v = 
at the redshift in Table 1). The error reported is the 
standard deviation of the Vcenter values found with the 
three fiux measurement methods and different channel 
averaging described in Section A. 2. The reported veloc- 
ity width (AV) is the full width of the emission, where 
'source' velocity channels are selected by eye. We give 
the velocity width of a single channel as the error. The 
line ratio rsi is given by 



rsi 



L' 



iL'c 



CO(3-2)/^CO(l-0) 



(3) 



and the error is calculated by propagating the errors of 
the individual line luminosities. 
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Name 


Transition 


(Jy km s"-*^) 


^CO 

(10^ K km s-i pc2) 


'^center 

(km s-i) 


Ay 

(km s-i) 


Non-matched r^-^ 


CI 


CO( J = 1^0) 
CO( J = 3^2) 


2.05±0.34 
9.02^±2.65 


8.25±2.83 
4.03"±1.69 


-84.9±18.1 
-49. 9± 7.3 


542.2±41.7 
542.3±23.4 


0.49«±0.26 


C2 


CO( J = 1^0) 
CO( J = 3^2) 


2.35±0.10 
8.39±1.12 


10.70±3.24 
4.25±1.40 


-6.7± 1.2 
19. 5± 8.2 


253.7±42.3 
169.2±42.3 


0.40±0.18 


C3 


CO( J = 1^0) 
CO( J = 3^2) 


4.39±0.09 
13.46±1.56 


21.70±6.53 
7.39±2.38 


-24.7± 5.7 
-14.0± 7.8 


384.0±42.3 
426.7±42.7 


0.34±0.15 


C4 


CO( J = 1^0) 
CO( J = 3^2) 


3.01±0.11 
11.88±1.22 


12.30±3.72 
5.41±1.72 


17.3±13.0 
33. 5± 5.0 


292.5±41.8 
250.7±41.8 


0.44±0.19 



Table 2 

Properties of the CO emission (described in Section 3.3). For each galaxy, the CO hne flux (Sco), luminosity (L^qq, see Equation 2), 
central velocity ('^center) and full velocity width (AV) are given for the CO(J = 1—^0) and CO(J = 3 — )► 2) lines. The ratio of the line 

luminosities (rsi) for each galaxy is given in the last column. 
^ CO( J = 3—^2) values for source CI are upper limits. 
^ rsi reported here is calculated from the total line luminosities (see Section 4.1). We perform a more careful calculation, matching the 1 

and 3mm data, in Section 4.3. 



Moment maps (discussed in Section 4.2 and displayed 
in Appendix B) are created from the 2a clip smooth mask 
(see Section A.2). Moment (total intensity) maps are a 
simple sum of the masked images in the 'source' velocity 
channels. Moment 1 (intensity- weighted mean velocity) 
maps are produced by summing the masked image mul- 
tiplied by the velocity in each channel then normalizing 
by the moment value. 

4. ANALYSIS 

We first present our total rsi measurements (using the 
total flux observed in each line) in Section 4.1. All four 
observed galaxies are discussed. However, the total rsi 
values use hne luminosities calculated completely inde- 
pendently of one another, making no attempt at match- 
ing the spectral or spatial resolution. In Section 4.2, we 
carefully match the C0( J =1^0) and C0( J = 3^2) 
data to each other as well as possible, in both spectral 
and spatial resolution. In Section 4.3 we re-calculate rsi 
using both the natural and adjusted spatial resolution, 
with matched spectral resolution in both cases. In Sec- 
tion 4.4 we investigate the radial dependence of rsi by 
fitting two-dimensional Gaussians to the total intensity 
maps in order to derive peak intensities, which are inde- 
pendent of spatial resolution. The analysis in Sections 
4.2-4.4 is performed only for the three galaxies that are 
detected significantly in both lines (CI is excluded due 
to marginal detection of the CO (J = 3 2) line). 

4.1. Total rsi 

As a first step, we calculate rsi in the simplest way 
possible, taking the ratio of the total line luminosities. 
Of the four galaxies observed, we detect three in both 
lines, yielding rsi values (see Table 2) of 0.34 to 0.44 
with an average value of 0.39. For source CI we derive 
an upper limit on rsi of 0.49, which is consistent with 
the other values. 

These rsi values carry the caveat that they are de- 
termined from line luminosities that are calculated in- 
dependently of one another. The 'source' velocity chan- 
nels for each line are selected by eye for the most con- 
vincing detection. Further, the CO (J =1^0) and 
CO (J = 3 2) lines are observed at different frequen- 
cies in the same configuration (except source C4). This 
yields different sampling of the ii^'-plane, which means 
that the two maps are sensitive to slightly different spa- 
tial scales. We re-calculate rsi more carefully for C2, C3 
and C4 in the following sections. 



4.2. Matching the 1 and 3 mm Data 

In this section, we attempt to standardize the fiux cal- 
culated for each line as much as possible by matching 
both spatial and spectral resolution. First, we use the 
same spectral resolution in both the CO (J = 1^0) and 
CO (J = 3^2) data, selecting the starting channel so 
that the image velocity channels line up exactly. 

Next, to investigate the effect of the different uv sam- 
pling of the two lines, we re-create the images using only 
data at uv distances present in both the 1 and 3 mm 
data (matched uv data): 8-44 kA, which corresponds to 
spatial scales of ~ 2.5''-13''. Further, we force all the im- 
ages to have a standard pixel size (0.99''for all except the 
CO (J = 3^2) line for C2 and C3 using all uv data, for 
which we use 0.33''pixels since these images have higher 
resolution). Figures 8 - 10 in Appendix B show images 
and spectra calculated within the standard source regions 
(discussed in Section A.2) using all uv data (top panel) 
and using matched uv data (bottom panel). 

This i^v-distance restriction makes the beam sizes more 
similar, but not necessarily the same due to different 
sampling within the allowed iiv-distance range. Table 
3 compares the measured fiuxes and beam sizes using 
all uv data and matched uv data. The change in the 
measured fiux is < 6% for sources C2 and C4, which is 
less than the expected measurement error in the total 
fiux. However, the change in fiux is more significant for 
source C3 in both lines. For the CO (J = 1^0) line, 
the matched iiv-coverage excludes the short uv spacings, 
which are sensitive to large-scale structure. Therefore, 
the presence of extended emission in source C3 would ex- 
plain the decrease in fiux in the matched uv data. This is 



Name, 
Trans. 


All uv data 


Matched uv data 


Flux 
% Diff. 




Beam 




Beam 


C2 1-0 
3-2 


2.30 
7.06 


4.90^^ X 4.08^^ 
1.70'' X 1.56'' 


2.15 
7.00 


4.54" X 3.75" 
3.20" X 2.84" 


-6.5% 
-0.8% 


C3 1-0 
3-2 


4.29 
13.75 


4.97" X 4.06" 
1.71" X 1.46" 


3.71 
11.61 


4.62" X 3.77" 
3.12" X 3.02" 


-13.5% 
-15.6% 


C4 1-0 
3-2 


2.95 
12.79 


4.99" X 4.02" 
3.33" X 2.41" 


2.76 
12.01 


4.62" X 3.83" 
3.55" X 3.08" 


-6.4% 
-6.1% 



Table 3 

Comparison of uv data selections. Measured CO line fluxes {Sqq 
in Jy km s"-*^) and beam sizes for images made using all uv data 

and matched uv data. The fluxes reported here are calculated 
using the 2cr smooth masking technique discussed in Section A.2. 
The percent difference in the measured flux is deflned as 

100(S'co, match - 'S'cO,all)/'S'cO,all- 
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1 



CO(l-O) all uv 
CO(l-O) mat uv 
CO(3-2) all uv 
CO(3-2) mat uv 



C2 




CO(l-O) all uv 
CO(l-O) mat uv 
CO(3-2) all uv 
CO(3-2) mat uv 



1 




1 



CO(l-O) all uv 
CO(l-O) mat uv 
CO(3-2) all uv 
CO(3-2) mat uv 



C4 



3^ 



-300 -200 -100 100 200 300 
velocity (km s"^ ) 



-300 



-200 -100 100 200 
velocity (km s"^ ) 



300 



-300 -200 -100 100 200 300 
velocity (km s"^ ) 



9 1.0 

J 0.8 

^ 0.6 

ro 

J 0.4 
II 

^ 0.2|- 
0.0. 



all uv 
mat uv 



o 1.0 



u 0.8 




300 -200 -100 100 200 300 
velocity (km s"^ ) 



0.0. 



all uv 
mat uv 



I _i 



-300 -200 -100 100 200 300 °-300 -200 -100 100 200 300 
velocity (km s"^ ) velocity (km s"^ ) 



Figure 4. L^qq and rsi profiles for the three detected sources. The top panels show the L^qq profiles for the CO (J = 1—^0) (red) and 
CO (J = 3 — )■ 2) (blue) transitions using all uv data (all uv, solid lines) and matched uv data (mat uv, dashed lines). The shading under 
each profile indicates which velocity channels are considered part of the source emission for each transition. The rsi profile is plotted in 
the bottom panel for both uv selections (all uv in solid black, matched uv in dashed magenta). The average rsi value for each profile is 
plotted as a horizontal black dotted (all uv data) or magenta dash-dotted (matched uv data) line. 

supported by the radial profile of the CO (J = 1^0) cu- 
mulative flux in Figure 6 in Appendix A, which shows the 
matched uv data agreeing with the all uv data at small 
radii and diverging at larger radii. The C0( J = 3 ^ 2) 
hue, on the other hand, looses long uv spacings (sen- 
sitive to small-scale structure) in the matched uv data. 
However, the iX'U-distance restriction excludes a sizeable 
fraction of the uv data, which decreases the signal to 
noise in the image so that less flux is recovered using 
the masking technique. A comparison of the unmasked 
CO (J = 3^2) line fluxes in the all uv and matched uv 
datasets for source C3 shows a drop of only 7.2%, which 
is within the expected measurement error. 

In the rest of our analysis, we consider the emission 
measured in both cases of uv coverage. 

4.3. Integrated rsi Velocity Profiles 

We now make a more careful measurement of rai, with 
the CO(J = 1^0) and CO(J = 3^2) flux measure- 
ments as well-matched to one another as possible. For 
each uv data selection (all uv^ matched uv), we calculate 
rsi in each channel where flux is measured in both lines. 

Figure 4 shows the integrated L^q and rsi velocity 
profiles for each of the three significantly detected galax- 
ies. Lqq as a function of velocity is plotted in the top 
panels for each transition (CO (J = 1 ^ 0) in red, 
CO (J = 3 ^ 2) in blue), for all uv data (solid lines) and 
matched uv data (dashed lines). Source velocity channels 
for each line (chosen by eye) are indicated by shading in 
the corresponding color. The bottom panels show the rsi 
profiles for all uv data (solid black line, with the average 
value indicated by the horizontal black dotted line) and 
the matched uv data (dashed magenta line, with average 



value shown by the horizontal magenta dash-dotted line). 
For galaxy C2, we have excluded the velocity channels at 
±86 km s~^ from this analysis. At —86 km s~^, we do 
not detect the CO (J = 3 2) line despite a 3a detec- 
tion in the CO(J = 1^0) line. At +86 km s-\ the 
CO (J = 3^2) image suggests a weak (2cr) detection at 
the center, but the measured flux is dominated by a 3cr 
noise spike close to this emission. Therefore, we use only 
the central three velocity channels in the CO (J = 3 2) 
transition. 

In galaxies C2 and C3, the rsi profiles appear well be- 
haved, remaining roughly constant in all velocity chan- 
nels. On the other hand, galaxy C4 shows significant 
deviations from a flat profile, with a very low rsi at ~ 
km s~^ and a very high value at ^ 40 km s~^. Due to the 
modest signal to noise ratio at which the CO (J = 3 2) 
line is detected in each channel, we cannot draw any 
definitive conclusions, but note that if the enhanced rsi 
is real, it may indicate a region of enhanced gas excita- 
tion (at a velocity close to 0, it is likely to be physically 
near the center of the galaxy), such as a starbursting 
clump, in which a higher rsi would be expected. 

For each galaxy. Table 4 gives the number of veloc- 
ity channels in which rsi is calculated, the mean and 
standard deviation (a) of rsi for each uv data selection 
and the the average rsi value for the two uv data selec- 
tions (galaxy mean) , with the expected measurement er- 
ror (meas. error). We expect a 20% measurement error 
in the flux of each line, in each channel (Section A. 2), 
which gives a 30% error in rsi in each channel, and a 
30%(A/'ch)~^^^ error in the average rsi- This value is re- 
ported in the last column of Table 4 for each galaxy. The 
bottom rows give the mean and standard deviation of the 
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Name 




All 


uv 


Matched uv 


Galaxy 


Meas. 






Mean 


a 


Mean 


a 


Mean 


Error 


C2 


3 


0.45 


0.06 


0.47 


0.14 


0.46 


0.08 


C3 


8 


0.39 


0.12 


0.39 


0.14 


0.39 


0.04 


C4 


6 


0.54 


0.32 


0.50 


0.32 


0.52 


0.06 


^^ample Mean 
Sample a 


U.46 
0.07 





Table 4 

Summary of EGNoG galaxy rsi values. For each source, we give 
the number of velocity channels in which rsi is calculated (Nch), 
the mean and standard deviation (a) of rsi using all uv data and 
the matched uv data, and the average rsi value for each galaxy, 
with the expected measurement error. In the bottom rows we 
give the mean and standard deviation of rsi in this sample of 
three galaxies. The expected measurement error for the average 
rsi in each galaxy is the expected measurement error in each 
channel (30%) divided by \/Nch- 



average rsi values for the three sample galaxies. In this 
paper, the standard deviation given is the square root of 
the unbiased sample variance. 

The average values for each uv selection are consistent 
with each other within each galaxy and roughly consis- 
tent across the three galaxies. While the uv selection ap- 
pears to have a small effect on the average rsi value, it is 
not in a systematic direction. The standard deviation of 
rsi we observe in galaxies C2 and C3 is consistent with 
the expected 30% measurement errors in each channel. 
The standard deviation for galaxy C4 is roughly twice 
the expected 30% error due to the two outlier channels 
discussed above. Overall, in the three EGNoG galaxies 
discussed here, we find a mean rsi of 0.46 with a stan- 
dard deviation of 0.07. Since we estimate each line flux 
has systematic errors of up to ^ 30% (see Section A. 2), 
the systematic errors in the ratio rsi will be < 40%. 

4.4. Radial Dependence of rsi 

Since rsi traces the local excitation conditions of the 
molecular gas of a galaxy, it is expected to vary within 
the disk of the galaxy. In fact, Dumke et al. (2001) found 
CO (J = 3 2) emission to be more centrally concen- 
trated than CO (J =1^0) emission in nearby galaxies, 
so that rsi decreased with radius. To look for radial vari- 
ation in rsi in the EGNoG data, we cannot compare the 
emission maps directly due to the marginal resolution 
of the galaxies and the different uv coverage of the two 
transitions. In order to disentangle the true emission dis- 
tribution from the uv sampling, we fit a two-dimensional 
Gaussian to the total intensity (moment 0) map of each 
CO transition using the MIRIAD program imf it. This 
analysis is presented in detail in Appendix C. 

While the error bars are large, we do find the ratio of 
the deconvolved peak intensities (representative of the 
conditions in the central region of the galaxy) to be sys- 
tematically higher than the ratio of the total fluxes. In 
our sample, we find rsi (peak) = 0.73 ±0.31, with a stan- 
dard deviation across the three galaxies of 0.04. We 
compare this to the ratio calculated from the fit total 
fluxes: 0.42 ± 0.11, with a standard deviation of 0.04. 
While the small standard deviations in our rsi values 
show consistency between the three galaxies, these re- 
sults are plagued by sizable uncertainties as a result of 
fitting data with only modest signal to noise (note that 
our error estimates are conservative and may in fact be 
over-estimating the true error). Therefore, we conclude 
(but not robustly) that rsi is higher in the center of the 



EGNoG galaxies than in the molecular disk as a whole 
(consistent with Dumke et al. 2001). 



5. DISCUSSION 



The observed rsi is a function of many parameters, 
including the excitation temperature (Tex), optical depth 
(r^y) and filling factor of each line (e.g. Hurt et al. 1993). 
For each line, the rest-frame intensity is 



--m (i - 

:5.(Tb) 



exp-"^) [5,(Te: 
_ 2/c^2^b,RJ 



^.(Tcmb)] (4) 



where Bjj{T) is the specific intensity of a blackbody at 
temperature T, Tb is the brightness temperature of the 
line, Tb,Rj is the Rayleigh- Jeans definition brightness 
temperature of the line and r^ is the efficiency of the 
coupling between the beam and source emission, which 
is determined by the beam size, the molecular gas disk 
size and the fraction of the disk that is emitting in the line 
(the filling factor). Note that Equation 4 is in the rest 
frame. The observed line intensity would be redshifted so 
that the corresponding observed brightness temperature 
(for both definitions) is reduced by a factor of (1 + z). 

For equal velocity widths in both lines (as we have en- 
forced in Sections 4.3 and 4.4), rsi, as defined in Equa- 
tions 2 and 3, reduces to the ratio of the Rayleigh- Jeans 
brightness temperatures. Following Equation 4, we can 
write rsi as 



rsi 



Tb,Rj(3-2) _ [^RJ^f^r^CMBTex] 



(3-2) 



Tb,Rj(l - 0) [/?Rj/?f/?r^CMBTex] 



(1-0) 



(5) 



where each term in Equation 4 is represented as a cor- 
rection factor ri so that we can express rsi more directly 
in terms of the line excitation temperatures. Specifi- 
cally, r^cMB gives the correction from Tex to the effec- 
tive excitation temperature (Tex,eff, where 5iy(Tex,eff) = 
Bu{Te^)-Bj,{TcMB)) due to the'CMB term. At z = 0.3, 
the temperature of the cosmic microwave background 
(Tcmb) is ~ 3.5 K, so that for Tgx =10-30 K, t/cmb ^ 
0.85 - 0.95 for the C0( J = 1^0) line and 0.98 - 1.00 for 
the CO (J = 3 2) line. The term accounting for the 
optical depth of each line is r]r = — e"'^^). The filling 
factor term, r^f (defined above) is expected to be < 1 and 
may be different for the two lines. Finally, we represent 
the difference between Tb and Tb rj by t^rj, which, for 
Tb = 10 - 30 K, is 0.75 - 0.91 for the C0( J =1^0) line 
and 0.39 - 0.75 for the CO(J = 3^2) line. Therefore, 
for Tex = 10 - 30 K for each line, we expect 



rsi = (0.4 -1.2) 



[^f^rT, 



exJs-2 



(6) 



exJi_o 



The excitation temperature of a given emission line 
will fall between the radiation temperature (Tcmb) and 
the kinetic temperature, depending on the density of the 
gas. Where the molecular gas density is larger than 
the critical density (for 20 K gas, ncrit ^ 10^ cm 

3 - 



CO(J = 1 ^ 0) and ^ 10^ cm"^ for CO(J 



3 for 
2)), 



the excitation temperature is equal to the kinetic tem- 
perature of the gas. Thus we see that optically thick 
gas with equal excitation temperature in both lines will 
have rsi « 1.0 ± 0.5 for Tex ^ 10 — 30 K (depending on 
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the relative filling factors and optical depths of the two 
lines). This is consistent with the general criterion that 
gas with rsi < 1 is sub-thermal and rsi > 1 is thermal- 
ized. Note that deviation from rsi = 1 in thermalized 
gas is a strong function of excitation temperature (due to 
^CMB and r]Rj): higher excitation temperature gas will 
have rsi closer to unity when thermalized. More gen- 
erally, if one considers various optical depths (assuming 
local thermodynamic equilibrium (LTE) with equal exci- 
tation temperatures), rsi greater than unity is indicative 
of warm, optically thin gas and a ratio less than unity is 
indicative of optically thick gas (e.g. Meier et al. 2001). 

Therefore, while we can comment on whether the exci- 
tation of the CO Jupper state is likely to be sub-thermal or 
thermalized, determining the conditions of the gas more 
accurately requires modeling of multiple transitions (e.g. 
using Large Velocity Gradient (LVG) models). Since one 
line ratio does not provide strong constraints, we do not 
perform LVG modeling as part of this work. In the fol- 
lowing sections we compare the rsi values found in the 
EGNoG galaxies to previous studies at low and high red- 
shift and we discuss the implications of this work for the 
interpretation of CO measurements in intermediate and 
high redshift galaxies. 

5.1. Comparison with Previous Work 

To date, the study of CO line ratios in intermediate and 
high redshift galaxies has been dominated by work on ex- 
treme starbursting systems: SMGs and quasars. The CO 
lines from ( J = 1 ^ 0) up to ( J = 9 ^ 8) in quasars at 
z ^ 2 — 4 are well-fit by a single component of highly ex- 
cited gas (Riechers et al. 2006; WeiB et al. 2007; Riechers 

2011) . In contrast, while the high- J CO transitions in 
z ~ 2 — 4 SMGs are fit by similar highly excited gas, re- 
cent observations of the CO (J = 1^0) line in these sys- 
tems reveal a diffuse, low-excitation component in addi- 
tion to the highly excited component (Carilli et al. 2010; 
Riechers et al. 2011), similar to what has been observed 
in local ULIRGS (Papadopoulos et al. 2007; Greve et al. 

2009) . In contrast, the limited work on z ~ 1 — 2 SFGs 
suggests low-excitation gas similar to the Milky Way 
(Dannerbauer et al. 2009; Aravena et al. 2010). Since 
the EGNoG galaxies are SFGs (lying on the main se- 
quence), we restrict the rest of this discussion to SFGs 
(normal star- forming galaxies) at low and high redshift. 

To place the rsi values we measure at z ~ 0.3 in 
the context of previous work on normal SFGs, Figure 5 
shows rsi versus approximate SFR for the EGNoG galax- 
ies (red diamonds) and a compilation of literature data, 
which includes three large surveys of rsi in nearby galax- 
ies (Yao et al. 2003; Mao et al. 2010; Papadopoulos et al. 

2012) as well as the study at z = 1.5 (Aravena et al. 

2010) . We plot rsi from our matched analysis (Section 
4.3) for sources C2, C3 and C4 (values from Table 4, with 
40% error bars to indicate potential systematic errors) 
and the upper hmit for source CI derived from the to- 
tal line luminosities (Section 4.1, Table 2). While SFRs 
are available for the present survey and Aravena et al. 
(2010), the other three surveys only provide total infrared 
luminosities (i^m) or far infrared luminosities (Lpm) cal- 
culated from IRAS (the Infrared Astronomical Satellite) 
fiuxes. We use approximate SFRs of 1.7 x 10~^^Lfir Mq 
yT~^ (i^FiR in L©; Kennicutt 1998) and LIR/Lfir = 1.3 
(Gracia-Carpio et al. 2008) in Figure 5. Values for the 
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Figure 5. Compilation of rsi literature data. rsi is plot- 
ted against approximate SFR (1.7 x 10~-'^*-'Lfir when a SFR 
is not available) for the following datasets: Yao et al. (2003) 
(Yao03), Mao et al. (2010) (MaolO), Papadopoulos et al. (2012) 
(Papal2), Aravena etal. (2010) (AravlO), this paper (EGNoG) 
and Fixsen et al. (1999) (MW GO, inner and outer). Milky Way 
points are shown at the left side of the plot for clarity. The black 
horizontal dashed line shows rsi = 1 and the horizontal solid lines 
indicate the average values for Yao03, MaolO and Papal 1 (the cor- 
responding shaded rectangles show the standard deviation around 
the average values). 

Milky Way are plotted as well: the galactic center (MW 
GC), inner disk (MWinner), and outer disk (MWouter) 
values, all in black, are taken from Fixsen et al. (1999) 
and plotted (slightly offset from one another) for illus- 
trative purposes at the left side of the plot. The aver- 
age rsi values of the Yao et al. (2003), Mao et al. (2010) 
and Papadopoulos et al. (2012) datasets are indicated by 
the horizontal lines of the corresponding color, with the 
corresponding shaded rectangles indicating the standard 
deviation around the average values. To guide the eye, 
the black dashed line shows rai = 1. 

Despite the large range of values observed, we find gen- 
eral agreement between our values, that of Aravena et al. 
(2010) and the average values of the local surveys, with 
the EGNoG galaxies showing slightly lower rsi values 
than the other studies. Note that the Mao et al. (2010) 
survey found an average value of 0.81 for the entire sam- 
ple, but a lower average of 0.61 for galaxies classified as 
'normal', based on SFR surface density as indicated by 
the far infrared luminosity and optical diameter. Fur- 
ther, the CO(J = 3^2) survey of local, low-SFR 
spirals by Wilson et al. (2012) find a lower value (using 
CO (J =1^0) luminosities from Kuno et al. 2007) of 
rsi ^ 0.2, but note that the authors have not made any 
correction for differences in the fraction of each galaxy 
mapped by the two surveys. Their rsi value is at the 
low end of the spread of values plotted in Figure 5. The 
range and average rsi values for the points plotted (in- 
cluding number of galaxies, redshift and SFR ranges) are 
given in Table 5. 
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Dataset 


N 


Redshift 
Range 


SFR^ 
(M0 yr-i) 


^31 


Range 


Mean 


Milky Way 


3 






0.4 - 0.9 


0.68 


Mao 2010 


61 


0.0 - 0.04 


0.03 - 200 


0.2 - 1.9 


0.81 


normal 


7 






0.3 - 1.5 


0.61 


starburst 


25 






0.3 - 1.9 


0.89 


Yao 2003 


60 


0.007 - 0.05 


0.1 - 100 


0.2 - 1.7 


0.66 


Papa 2012 


70 


0.006 - 0.08 


0.7 - 150 


0.1 - 1.9 


0.67 


EC^NoC^ 


3 


0.28 - 0.31 


39 - 65 


0.39 - 0.52 


0.46 


Arav 2010 


1 


1.52 


220 




0.61 



Table 5 

Summary of rsi from the literature. The range and mean values 
are given for the following datasets: Fixsen et al. (1999) (Milky 
Way), Mao et al. (2010) (Mao 2010), Yao et al. (2003) (Yao 
2003), Papadopoulos et al. (2012) (Papa 2012), this paper 
(EGNoG) and Aravena et al. (2010) (Arav 2010). Columns 2 to 4 

give the number of galaxies, redshift and SFR ranges for each 
dataset. For Mao 2010, we also report the rsi range and average 
for the 'normal' and 'starburst' subsets. 
^ SFR values are approximate (see Section 5.1). 



In order to compare our values to the literature data, 
we first discuss the nature of the observations presented. 
Yao et al. (2003) observe 60 IR-luminous galaxies (most 
sources have Lfir > 10^^ L©) selected from the SCUBA 
Local Universe Galaxy Survey (SLUGS). Each line is 
measured with a single pointing of a single-dish telescope: 
CO(J = 1 0) at the Nobeyama Radio Observatory 
(NRO), CO(J = 3 ^ 2) at the James Clerk Maxweh 
Telescope (JCMT). For both measurements, the beam 
size is ~ 15'', corresponding to physical sizes of ~ 0.5 — 13 
kpc for their sample galaxies (~ 0.5 — 5 kpc for most of 
the sample). 

Mao et al. (2010) measure the CO(J = 3 ^ 2) line for 
125 nearby galaxies of various types (e.g. normal, star- 
burst, LIRG, ULIRG) using the Heinrich Hertz Telescope 
(HHT; beam size ^ 22''). This beam size corresponds to 
a physical size of ~ 0.25 — 17 kpc for the sample galaxies 
(1.7 kpc on average). The CO(J = 1^0) data for their 
sample were taken from the literature and therefore the 
beam size depends on the telescope used. In their anal- 
ysis, Mao et al. (2010) only use those galaxies for which 
they found IRAM 30m CO(J = 1^0) data (61), which 
have a beam size matching the HHT CO (J = 3 2) 
data. These measurements are plotted with error bars in 
Figure 5. Galaxies with C0( J = 1^0) data from other 
sources are reported as upper or lower limits. 

Papadopoulos et al. (2012) examine a composite sam- 
ple of 70 LIRGs in the nearby universe {z < 0.1) spanning 
a wide range of morphologies. They present new mea- 
surements of 36 galaxies and data from the literature for 
34 more. The new measurements use the IRAM 30m 
telescope for CO (J =1^0) (beam size ~ 22") and 
the JCMT for CO(J = 3^2) (beam size ^ 14"), and 
therefore are not observed with the same beam size. The 
authors do not comment on matching beam sizes in the 
new measurements or those from the literature. 

The galaxy in which Aravena et al. (2010) measure rai, 
BzK-21000 at ^ = 1.5, was observed in the CO(J = 
3^2) line with the Plateau de Bure Interferometer 
(PdBI) by Dannerbauer et al. (2009) and in the C0( J = 
1^0) line with the Very Large Array (VLA) in C and 
D configurations (Aravena et al. 2010). In both cases, 
the source is unresolved or marginally resolved, so the 
measured fluxes should represent the total emission of 
the galaxy. 



The three samples of nearby galaxies are observed 
with single-dish instruments, typically sampling the in- 
ner portion of the molecular gas disk of the observed 
galaxies. This is different from the rsi measurements 
for the EGNoG galaxies (Table 4) and the galaxy at 
z = 1.5 (Aravena et al. 2010), for which the emission 
from the entire gas disk is olDserved. We note that a 
gradient in rsi has been observed in nearby galaxies 
(Dumke et al. 2001), with higher values measured in the 
centers. This can be inferred from the Milky Way mea- 
surements (Fixsen et al. 1999) as well. This effect may 
well account for the EGNoG and Aravena et al. (2010) 
measurements appearing systematically lower than the 
average values reported by the surveys observing the cen- 
tral regions of nearby galaxies. Supporting this explana- 
tion, our analysis in Section 4.4 suggests rsi is higher 
in the centers of the EGNoG galaxies (more in line with 
the average values of the local surveys) compared to rsi 
averaged over the entire disk. However, with such a large 
spread in values and inhomogeneous datasets, it is diffi- 
cult to draw any robust conclusions on this point. 

5.2. Implications 

We have measured rsi = 0.46 ±0.07 (with systematic 
errors of up to 40%) in three galaxies at 2; « 0.3, sug- 
gestive of optically thick, sub-thermal gas. Despite being 
massive and highly star-forming (with SFRs of 50-65 M© 
yr~^ and stellar masses of ^ 2 x 10^^ M©), the excita- 
tion of the gas in these galaxies is consistent with SFGs 
like local spirals, not star-bursting systems like ULIRGS, 
SMGs and quasars. Since the EGNoG galaxies have been 
selected from the main sequence of star-forming galaxies 
at z = 0.3, our flndings suggest that galaxies on the main 
sequence, over a range of star formation activity, harbor 
sub- thermally excited gas. Therefore, we suggest that 
CO line ratios similar to those observed in local spiral 
galaxies are appropriate for main-sequence star-forming 
galaxies. 

The extension of the EGNoG results to main-sequence 
star-forming galaxies means that sub-thermal line ra- 
tios are appropriate for the z ^ 1 — 2 SFGs in which 
Tacconi et al. (2010) and Daddi et al. (2010) report high 
molecular gas fractions (20-80%). Tacconi et al. and 
Daddi et al. estimated the molecular gas mass associ- 
ated with the observed Jupper > 1 CO luminosity assum- 
ing sub-thermal rji line ratios and a Milky Way-like aco 
(Daddi et al. 2010 used a slightly smaller value for aco)- 
As the conversion factor is a complex problem, expected 
to vary as a function of gas excitation, density, metal- 
licity and radiation fleld (Shetty et al. 2011; Leroy et al. 
2011), the results of this work do not directly inform 
the choice of conversion factor. However, the rsi line 
ratios of the EGNoG gas excitation sample support the 
assumption of a sub-thermal line ratio in these studies. 
Speciflcally, Tacconi et al. (2010) use rsi = 0.5 to calcu- 
late molecular gas masses, which is consistent with our 
results (Daddi et al. 2010 observed the CO(J = 2^1) 
line and assume a sub- thermal r2i). 

6. CONCLUSIONS 

This paper presents the gas excitation sample of the 
EGNoG survey. We report robust detections of the 
C0( J = 3^2) and C0( J =1^0) lines in three galax- 
ies at z ^ 0.3, and an upper limit for the fourth galaxy. 
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The average rsi value for this sample is 0.46 ±0.07, with 
systematic errors of up to 40%. This value is consistent 
with published r^i values for main-sequence star-forming 
galaxies at 2; ~ as well as the single measurement of at 
z > 0.3 (Aravena et al. 2010). The sub-thermal excita- 
tion of the CO (J = 3 2) hue suggests the excitation 
state of the molecular gas in these galaxies is similar to 
local spirals, and is not indicative of a starburst. We 
conclude that the galaxies in our sample (and by exten- 
sion, the main sequence galaxies at z 1 — 2 studied by 
Tacconi et al. 2010 and Daddi et al. 2010) harbor cold, 
optically thick molecular gas despite being massive and 
highly star-forming. 
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APPENDIX 

A. DATA REDUCTION AND FLUX MEASUREMENT 
A.l. Data Reduction 

Each dataset is reduced and calibrated as follows. The data are flagged for antenna - antenna shadowing as well 
as any other issues during the observation. The instrument bandpass is calibrated with mf cal on a bright passband 
calibrator. The time-dependent antenna gains (from atmospheric variation) are derived by performing a self cal on 
the phase calibrator with an averaging interval of 18 minutes (the timescale of switching between the source and phase 
calibrator). For sources CI and C2, the phase calibrator, 0854+201, is ~ 10% polarized, which required additional 
steps in the reduction of the 3mm data (observed with linearly polarized feeds). This is described in more detail in 
the full survey paper, Bauermeister et al. 2013b, in preparation. 

For each dataset, the flux of the phase calibrator is set during the antenna gain calibration in order to properly 
set the flux scale of the data. The flux of each phase calibrator is assumed to be constant over timescales of weeks, 
and is therefore determined from the best datasets of the survey using boot flux on bandpass-calibrated, phase-only 
gain-calibrated data with Mars or MWC349 as a primary flux calibrator. The brightness temperature of Mars is set 
by the CARMA system using the Caltech thermal model of Mars (courtesy of Mark Gurwell), which includes seasonal 
variations in temperature and can be accessed in MIRIAD using marstb. This model gives brightness temperatures 
^ 218 K at 266 GHz and ^ 208 - 200 K at 88 GHz for August-November 2011 observations, and ^ 201 K at 266 GHz 
for April 2012 observations. For MWC349 we set the flux to 1.2 Jy at 88 GHz, the typical value from historical flux 
monitoring at CARMA. The fluxes used for each phase calibrator are as follows: for the 3mm data (^ 89 GHz, August 
to November 2010), the flux of 0854+201 is set to 4 Jy (9% linearly polarized), 1357+193 0.8 Jy and 1310+323 1.7 
Jy. For the 1mm data (^ 265 GHz, August 2011 and April 2012) the flux of 0854+201 is set to 2 Jy in August 2011 
and 4 Jy in April 2012, 1224+213 0.6 Jy in April 2012 and 1310+323 0.6 Jy in August 2011. 

Images are produced combining all fully calibrated datasets for each source. We used invert, weighting the visibilities 
by the system temperature as well as using a Briggs' visibility weighting robustness parameter (Briggs 1995) of 0.5. 
Since CARMA is an inhomogeneous array (these data use both the 10m and 6m dishes), we also use options=inosaic 
in the invert step in order to properly handle the three different primary beams patterns (lOm-lOm, 10m-6m and 
6m-6m). All observations are single-pointing. The resulting images are primary-beam-corrected. In most of the sources 
discussed here, we map the 3mm data (CO (J =1^0) line) channel by channel, using the full spectral resolution of 
42 km s~^. We match this resolution in the 1mm data (CO(J = 3^2) line) by averaging channels in sets of three. 
The exception to this scheme is source CI, which has emission spread over a very wide velocity range, requiring more 
channel averaging to increase the signal to noise. 

We deconvolve each image with mossdi (the mosaic version of clean), cleaning down to the rms noise within a single 
channel, within a cleaning box selected by eye to include only source emission. We clean only channels which contain 
visible source emission. Cleaning down to a specifled noise level is preferred to using a set number of clean iterations 
due to the nature of the spectral line emission: some channels will contain more flux and therefore require more clean 
iterations. In tests using a model source of known flux inserted into real data (emission- free channels), we found a la 
cutoff to best extract the true source emission without overestimating the flux over a range of detection signiflcance 
levels, with a 10-30% uncertainty in the recovered flux (depending on the significance of the signal). The final clean 
images are produced by rest or, which convolves the clean component image with the clean beam (calculated by fitting 
a Gaussian to the combined mosaic beam given by mospsf ), and adds the residuals from the cleaning process. 

A. 2. Flux Estimation 

Total source fiuxes in the CO lines are calculated by summing 'source' pixels in each 'source' velocity plane of the 
image. The source velocity planes are selected by eye. The source pixels are those within the 'source' region which 
are not masked by our smooth mask (adapted from Dame 2011). The smooth mask is created by applying a 2cr clip 
to a smoothed version of the image: Hanning smoothing is done along the velocity axis, and each velocity plane is 
convolved with a Gaussian beam twice the size of the original synthesized beam. This smoothed mask is used in order 
to exclude noise pixels but still capture low- level emission that a simple clip would miss. In our own testing (using a 
model source of known fiux inserted into real data), we found a 2a clip to best reproduce the true fiux over a range of 
detection significance levels. 

The appropriate source region size is selected to recover all of the fiux without including the negative bowl. Since 
we do not have single dish data to complement the interferometric data presented here, our datasets are missing the 
shortest uv spacings. As a result, the emission in the clean images sits in a negative bowl, which will affect the 
measured fiux (calculated by summing pixels within a given radius). In order to accurately estimate the fiux, we use 
the radius at which the radial profile of the enclosed fiux first peaks, thereby excluding the negative bowl. The radial 
profiles of the enclosed fiux for sources C2, C3 and C4 are shown in Figure 6: C0( J = 1 ^ 0) in red in the top panels 
and CO(J = 3 ^ 2) in blue in the bottom panels. The profiles for both uv selections (all and matched), unmasked 
and masked are shown for each transition. The negative bowl is most evident in the unmasked profiles (dotted and 
dash-dotted lines), in which the enclosed fiux peaks and then decreases with radius. The radius of the first peak of 
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radius (") radius (") radius (") 

Figure 6. Enclosed flux (.Sco in Jy km s"-*^) as a function of radius for sources C2, C3 and C4. The top panels show the CO(J = 1 — )► 0) 
line flux for four cases: all uv data, masked (solid) and unmasked (dotted); matched uv data, masked (dashed) and unmasked (dash-dotted). 
The bottom panels show the CO (J = 3—^2) line flux for the same four cases. The vertical black dotted lines indicate the radius at which 
the enclosed flux peaks (6.5'' for CO(J = 1^0), 4.5" for CO(J = 3^2)). 

the enclosed flux distribution is 6.5'' for the CO (J = 1^0) data and 4.5'' for the CO (J = 3^2) data (shown by 
the vertical black dotted lines). These radii are used in the calculation of total fluxes throughout this work. 

The error in the flux measurement is estimated from the standard deviation of the measured fluxes using different 
velocity channel averaging and starting channel, using three different methods of calculating the flux in each case. The 
three methods are: the 2(j masking technique described above, the same masking technique with a 3cr clip, and the 
simple addition of all pixels (no mask) within the source region. 

We performed extensive testing of our analysis technique in order to choose the parameters of the reduction to 
eliminate systematic offsets and minimize the uncertainty due to noise (as described above). We find a 10-30% error 
in the flux measurement coming from noise in the data reduction and analysis steps. From this, we take an average 
uncertainty of 20% in the flux estimated in each channel, which results in a uncertainty in the total flux of 20%(A^ch)~^"^ 
{Nch is the number of velocity channels in which the flux is summed). For the total flux values reported in Table 2, 
we estimate the error from the variation in the flux calculated using different channel averaging, flux measurement 
method, etc. (described above), which is consistent with the 20%(A^ch)~^"^ we expect. In our analysis of the integrated 
flux velocity profiles (Section 4.3), we assume errors of 20% in the flux in each channel. 

Further, these data suffer from systematic errors due to absolute flux calibration and primary beam correction. We 
set the flux scale in our dataset based on a primary flux calibrator (Mars or MWC349), the flux of which is only known 
to ^ 20%. In the primary beam correction of the dataset, pointing and focus errors at the time of the observations 
as well as errors in the primary beam model can significantly reduce image fidelity, leading to errors in the measured 
fluxes of ^ 20% (see SKA Memo #103, Wright & Corder 2008). Combining these systematic errors in quadrature, 
we estimate that our flux measurements suffer from systematic uncertainties of up to ^ 30%. We consider all these 
factors in the presentation of our data in Table 2: for the line flux (^co), we report the measured error; for Lqq, we 
include a 30% systematic error (added in quadrature to the measured error in Sco)- 
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Figure 7. CO(J = 1 — 0) emission in source CI. The left panel shows the optical image. The moment and moment 1 maps are 
displayed in the left middle and right middle panels respectively. The dotted white ellipse indicates the source region (6. 5^^ radius). In the 
moment map, the beam size is indicated by the solid white ellipse in the lower left corner and a 10 kpc scale bar is given in the top right. 
The far right panel shows the spectrum of the galaxy: the solid blue line is calculated with masking, the dotted green line is without. 



B. MOMENT MAPS 

In this appendix, we present the moment maps, optical images and spectra of the CO emission in the EGNoG bin 
C galaxies. Figure 7 shows the CO(J =1^0) emission in galaxy CI (we do not detect C0( J = 3^2) emission in 
galaxy CI). Figures 8-10 show both CO(J = 1^0) and CO(J = 3^2) emission in galaxies C2, C3 and C4, using 
all uv data and matched uv data. 

C. RADIAL DEPENDENCE OF i^si 

Since rai traces the local excitation conditions of the molecular gas of a galaxy, it is expected to vary within the disk 
of the galaxy. In fact, Dumke et al. (2001) found CO(J = 3^2) emission to be more centrally concentrated than 
CO (J = 1^0) emission in nearby galaxies, so that rsi decreased with radius. To look for radial variation in r^i in 
the EGNoG data, we cannot compare the emission maps directly due to the marginal resolution of the galaxies and the 
different uv coverage of the two transitions. We emphasize that the uv coverage (and thus different spatial resolution) 
of each transition determines the shape of the radial profile of the enclosed flux (see Figure 6 in Appendix A), making 
a direct ratio of the two radial profiles meaningless without perfectly matched uv sampling (even our matched uv data 
do not meet this criterion due to different sampling within the allowed iX'U-distance range). 

In order to disentangle the true emission distribution from the uv sampling, we fit a two-dimensional Gaussian to 
the total intensity (moment 0) map of each CO transition using the MIRIAD program imf it. The program fits for the 
position, peak intensity, total flux, size, and deconvolved source size. We perform the Gaussian fit on four versions of 
the moment map for each transition in each galaxy: for each uv data selection (all and matched), we use moment 
maps produced with and without 2cr smooth masking. In order to derive a good fit to the source emission, we restrict 
the Gaussian fit to the standard source region (Section A. 2). The error in the fit total flux is calculated by propagating 
the errors on the fit peak intensity and (not deconvolved) size. Since these parameters are not independent, this error 
is likely an overestimate of the true error in the fit total flux. 

The residual of each fit is inspected. We note that in general, the emission is fairly well fit by a Gaussian except for 
the CO (J = 3^2) transition in galaxies C2 and C3, which are observed at higher resolution. The structure present 
it in these higher resolution images is not well-fit by a Gaussian. 

We look for a radial dependence in r^i by comparing the ratio calculated from the fit total fluxes to the ratio 
calculated from the deconvolved peak intensities. The results are presented in Table 6. For each galaxy, for each 
combination of uv data selection (all, matched) and masking (no mask, masked), we derive the deconvolved peak 
brightness temperature (Tb) for each CO transition and calculate rai from the deconvolved peak intensities as well as 
the fit total fluxes (all values and errors are derived from the parameters of the Gaussian fitting) The rest-frame peak 
brightness temperature (in K) is given by 

n = . (CD 



/cb In 



2hv 



c2(l + ^)3/p. 



where is the rest frequency of the transition and /peak is the observed peak specific intensity (erg s ^ cm ^ Hz 
Sr~^). We calculate /peak from the velocity width (AV in km s~^), fit total flux {Sqq in Jy km s~^) and fit deconvolved 
size (FWHMmajor and FWHMminor in arcseconds): 

^ 4(ln2)(206265^)5co .^2) 

peak ^(1023)AyFWHM^ajorFWHM„,inor ^ ^ 

The ratio of the lines, rai, is calculated from the line flux (total or peak) according to Equations 2 and 3. Note that 
this ratio is not equivalent to the ratio of the brightness temperatures as defined here (see discussion in Section 5). 

The errors in the quantities given in Table 6 are calculated from the errors in the relevant fit parameters reported by 
imf it. We estimate the error in the deconvolved source sizes (FWHMmajor and FWHMminor) from the spread in the 
four estimates (all and matched uv^ with and without masking) for each transition, for each galaxy. For each quantity. 
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Figure 8. Two sets of 6 panels showing the optical image (top left), moment maps (middle and right panels) and spectra (bottom left) 
for source C2: the top set uses all uv data, the bottom set uses matched uv data. In each set, the top middle and right panels show the 
moment (total intensity) and moment 1 (intensity- weighted mean velocity) maps, respectively, for the CO (J = 1—^0) line. The bottom 
middle and right panels show the same, but for the CO (J = 3 — )► 2) line. The synthesized beam (solid white ellipse) and 10 kpc scale 
(white bar) are indicated in each middle panel. The dotted white circles in the middle and right panels show the source regions in which 
flux is summed (with radii of 6. 5^^ and 4.5^^for the CO( J = 1 — )► 0) and CO( J = 3—^2) data respectively; see Section A. 2). The bottom left 
panels show the spectra in mJy, with the CO(J = 3—^2) spectrum reduced by a factor of 4.5 to match the scale of the CO(J = 1—^0) 
spectrum. 
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Figure 9. Same as Figure 8, but for source C3. 
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Figure 10. Same as Figure 8, but for source C4. 
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we report the average value for each galaxy, taking the error to be the average fractional error. From the average 
values, we report the mean (with an uncertainty from the errors in the average values) and standard deviation for the 
sample as a whole. 

We find peak brightness temperatures of approximately 2 K and 4 K for CO (J = 1^0) and CO (J = 3^2) 
respectively. This is lower than the excitation temperature observed in Giant Molecular Clouds (GMCs) in the Milky 
Way (e.g. 10-30 K, Polychroni et al. 2012), as expected since the filling factor for GMCs in the molecular gas disk is 
less than unity. 

While the value of rsi is a complex problem (discussed in detail in Section 5), the larger brightness temperature 
of the CO (J = 3^2) line (relative to the CO (J = 1^0) line) may be explained simply by differing excitation 
temperatures or filling factors. For instance, if the CO (J = 1^0) traces an additional diffuse, lower-excitation gas 
component not traced by the CO (J = 3^2) line (as observed in SMGs by Carilli et al. 2010; Riechers et al. 2011), 
the average brightness temperature of the CO (J = 3^2) line will be higher than the CO (J = 1^0) line. On the 
other hand, a larger brightness temperature in the CO (J = 3^2) line may be due to a larger filling factor for the 
CO (J = 3^2) emission. Our Gaussian fits find systematically smaller sizes for the CO (J = 3^2) emission relative 
to the CO (J = 1^0) sizes, which is consistent with the results of Dumke et al. (2001), who found CO (J = 3^2) 
emission to be more centrally concentrated than CO (J = 1^0) emission in nearby galaxies. Since the molecular 
gas disks of galaxies tend to have an exponential radial profile (e.g. Regan et al. 2001; Leroy et al. 2009), the gas-rich 
central region of the galaxy will account for a larger fraction of the CO (J = 3^2) emission area, increasing the 
effective filling factor of the CO (J = 3^2) emission relative to that of the CO (J = 1^0) emission. Assuming 
a constant excitation temperature for both transitions, the larger effective filling factor of CO(J = 3 2) emission 
would produce a higher brightness temperature, as we observe. 

While the error bars are large, we do find the ratio of the deconvolved peak intensities (representative of the conditions 
in the central region of the galaxy) to be systematically higher than the ratio of the total fiuxes. In our sample, we 
find r3i(peak) = 0.73 ± 0.31, with a standard deviation across the three galaxies of 0.04. We compare this to the 
ratio calculated from total fiuxes: 0.42 ±0.11, with a standard deviation of 0.04. While the small standard deviations 
in our rsi values show consistency between the three galaxies, these results are plagued by sizable uncertainties as a 
result of fitting data with only modest signal to noise (note that our error estimates are conservative and may in fact 
be over-estimating the true error). Therefore, we conclude (but not robustly) that rsi is higher in the center of the 
EGNoG galaxies than in the molecular disk as a whole (consistent with Dumke et al. 2001). 



Na- 


uv 


Gaussian Fit Peak 


Gaussian Fit rsi 


me 


Sel. 


CO(l-O) 


CO(3-2) 


Peak 


Total 


C2 


all 


1.83 ±0.32 


3.79 ±0.40 


0.75 ±0.54 


0.41 ±0.26 






1.83 ±0.32 


3.71 ±0.34 


0.68 ±0.47 


0.48 ±0.29 




mat 


1.91 ±0.13 






0.37 ±0.13 






1.90 ±0.15 


3.86 ±0.47 


0.72 ±0.42 


0.42 ±0.20 




avg 


1.87 ±0.23 


3.79 ±0.40 


0.72 ±0.48 


0.42 ±0.22 


C3 


all 


1.78 ±0.32 


3.45 ±0.35 


0.53 ±0.40 


0.44 ±0.14 






1.65 ±0.27 


3.51 ±0.51 


0.74 ±0.67 


0.38 ±0.22 




mat 


1.68 ±0.28 


3.75 ±0.37 


0.94 ±0.67 


0.35 ±0.07 






1.65 ±0.27 


3.65 ±0.38 


0.88 ±0.66 


0.36 ±0.11 




avg 


1.69 ±0.28 


3.59 ±0.40 


0.77 ±0.60 


0.38 ±0.13 


C4 


all 


2.03 ±0.60 


4.10 ±0.40 


0.76 ± 0.72 


0.48 ±0.38 






2.09 ± 0.33 


4.16 ±0.44 


0.74 ±0.46 


0.45 ±0.16 




mat 


2.08 ± 0.36 


3.89 ±0.31 


0.57 ±0.34 


0.50 ±0.15 






2.23 ±0.37 






0.45 ±0.19 




avg 


2.11 ±0.42 


4.05 ±0.38 


0.69 ±0.50 


0.47 ±0.22 




mean 


1.89 ±0.19 


3.81 ±0.23 


0.73 ±0.31 


0.42 ±0.11 




a 


0.21 


0.23 


0.04 


0.04 



Table 6 

Gaussian fit Tb and rsi. Deconvolved, rest-frame peak brightness temperatures (peak Tb) and rsi calculated from the deconvolved peak 
intensities (rsi peak) and from the total fiuxes (rsi total). All values are calculated using the parameters obtained from the Gaussian 
fitting. For each galaxy, we report five values in each column: one for each of the four data selections used (all uv (all) and matched uv 

(mat), without masking first and with masking second) and the average value with a typical error. From the average values for each 
galaxy, we calculate the mean (with an uncertainty from the errors in the average values) and the standard deviation for each quantity, 

presented in the bottom rows. 



